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An  analytical  study  of  the  effect  of  diffusioosmosis  caused  by  the  concentration  gradient  of  hydrogen 
ions  on  the  isothermal  transport  of  water  in  a  fully  hydrated  membrane  of  a  polymer  electrolyte  fuel 
cell  (PEFC)  is  presented.  A  capillary  tube  or  slit  with  a  negatively  charged  wall  is  chosen  to  model  the 
nanopores  of  the  membrane.  The  electric  double  layer  adjacent  to  the  capillary  wall  may  have  an  arbitrary 
thickness  relative  to  the  capillary  radius  and  its  electrostatic  potential  distribution  is  determined  as  the 
solution  of  the  Poisson-Boltzmann  equation.  Solving  a  modified  Navier-Stokes  equation,  the  fluid  velocity 
in  the  axial  direction  of  the  capillary  induced  by  the  macroscopic  electric  field  and  protonic  concentration 
gradient  is  obtained  as  a  function  of  the  radial  position  in  closed  forms.  The  results  for  the  local  and 
averaged  electrokinetic  velocities  in  the  capillary  show  that  the  effect  of  diffusioosmosis  on  the  water 
transport  in  the  membrane  of  a  PEFC  can  be  significant  in  comparison  with  that  of  electroosmosis  under 
low-potential-difference  operations. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  polymer  electrolyte  fuel  cell  (PEFC)  uses  a  polymeric  mem¬ 
brane  as  an  electrolyte,  analogous  in  acidity  to  the  electrolyte  in 
the  automotive  battery,  but  dimensionally  fixed.  This  facilitates 
sealing  in  the  production  process  and  furnishes  both  cell  and 
stack  with  longevity.  PEFCs  can  generate  high  electric-current  den¬ 
sity  and  allow  quick  start-ups  and  instant  response  to  changes 
in  the  demand  of  power.  Being  efficient,  environmentally  clean, 
and  potentially  low-cost  power  systems  for  both  stationary  and 
mobile  applications,  PEFCs  are  of  interest  to  both  electric  vehicle 
and  portable  electronics  manufacturers  [1-3]. 

The  polymer  electrolyte  membrane  (PEM),  which  prevents  mix¬ 
ing  of  feed  gases  and  provides  transport  of  protons  from  the  anode 
to  the  cathode,  is  one  of  the  key  components  in  PEFC  systems.  Per- 
fluorosulfonic  acid  (PFSA)  copolymer  membranes,  such  as  DuPont’s 
Nation  (composed  of  carbon-fluorine  backbone  chains  with  perflu- 
oro  side  chains  terminated  by  negatively  charged  sulfonic  groups), 
are  the  most  commonly  used  PEMs  in  PEFCs  operating  at  low 
temperatures  below  80  °C.  Some  high-temperature  PEMs  for  opera¬ 
tion  above  100  °C,  such  as  naphthalene-based  sulfonated  polyimide 
(SPI)  and  modified  PFSA  membranes,  are  also  under  active  devel¬ 
opment  [4,5]. 
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The  performance  of  a  PEFC  is  strongly  influenced  by  the  protonic 
conductivity  of  its  PEM,  and  the  conductivity  is  highly  dependent 
on  the  state  of  hydration  of  the  membrane.  The  sources  of  water  in 
the  membrane  are  the  humidified  reactant  gases  and  the  product  of 
the  chemical  reaction  of  hydrogen  ions  and  oxygen  occurred  at  the 
cathode.  Undoubtedly,  the  protonic  conductivity  of  a  PEM  increases 
with  its  content  of  water.  But,  an  excess  of  water  in  the  fuel  cell 
may  lead  to  electrode  flooding  and  thus  the  diffusion  of  the  reactant 
gases  will  be  suppressed.  In  addition,  external  humidification  of  the 
feed  gases  is  a  burden  for  the  fuel  cell.  Therefore,  understanding  of 
the  water  transport  and  successful  management  of  the  water  in  the 
PEM  is  critical  to  optimize  the  performance  of  PEFCs  [6-8]. 

The  water  distribution  throughout  a  PEFC  is  determined  by  the 
interplay  among  the  water  uptake  by  the  PEM,  the  rate  of  chemical 
reactions  at  the  electrodes,  the  diffusion  of  water  due  to  its  activity 
gradients,  and  the  water  transport  in  the  membrane  accompanied 
by  the  protonic  current  from  the  anode  to  the  cathode.  The  protonic 
current  produces  water  flow  along  the  charged  nanopore  walls  in 
the  membrane  by  two  electrokinetic  mechanisms.  The  first  is  elec¬ 
troosmosis  resulting  from  the  interaction  between  the  macroscopic 
electric  field  and  the  electric  double  layer  adjacent  to  the  pore  wall, 
and  the  effects  on  the  water  transport  caused  by  this  well-known 
mechanism  were  studied  extensively  in  the  past  [9-14].  The  second 
mechanism  is  diffusioosmosis  (also  known  as  capillary  osmosis) 
comprising  the  stresses  developed  by  the  tangential  gradient  of  the 
excess  pressure  within  the  double  layer  [15-21  ],  and  this  effect  on 
the  water  transport  in  the  PEM  has  not  been  investigated  yet. 


712 


H.J.  Keh,  H.C.  Ma  /  Journal  of  Power  Sources  180  (2008)  711-718 


Nomenclature 

e  elementary  charge  (C) 

E  macroscopic  electric  field  (V  m_1 ) 

E  =E  e/ickT 

h  half  thickness  of  a  slit  (m) 

k  Boltzmann’s  constant  (J  I<-1 ) 

L  length  of  a  capillary  (m) 

n  concentration  distribution  of  hydrogen  ions  (m-3 ) 

n°°  bulk  concentration  distribution  of  hydrogen  ions 

(m-3) 

N  =|vnoo|//rnoo(0) 

p  pressure  distribution  (N  m-2 ) 

r  radial  coordinate  in  a  tube  (m) 

R  radius  of  a  tube  (m) 

T  absolute  temperature  (I<) 

u  fluid  velocity  distribution  (m  s-1 ) 

(u)  average  fluid  velocity  (m  s-1 ) 

if  characteristic  electrokinetic  velocity  (m  s-1 ) 
y  lateral  coordinate  in  a  slit  (m) 

z  axial  coordinate  in  a  capillary  (m) 

Greek  letters 
y  =ln(/cr) 

s  dielectric  permittivity  (C2  J-1  m-1 ) 

f  zeta  potential  of  capillary  walls  (V) 

i  -teller 

rj  fluid  viscosity  (kg  m-1  s-1 ) 

k  Debye  screening  parameter  (nrr1 ) 

v  =  —xjr  +  2  ln(/cr) 

o  surface  charge  density  (C  m-2 ) 

x/r  electrostatic  potential  distribution  (V) 

\j/  =_\lrelkT 

xj/0  \j/  at  the  axis  of  a  capillary 


The  fluid  flow  resulting  from  the  interaction  between  a  solid  wall 
and  a  tangential  solute  concentration  gradient  is  termed  diffusioos- 
mosis.  The  ionic  solute-charged  wall  interaction  in  diffusioosmosis 
is  electrostatic  in  nature  and  its  range  is  the  Debye  screening  length 
/c-1  (defined  by  Eq.  (4)).  The  fluid  motion  caused  by  diffusioosmo¬ 
sis  has  been  analytically  examined  for  ionic  solutions  near  a  plane 
wall  [22,23]  and  inside  a  capillary  pore  [23-27].  Some  experimen¬ 
tal  results  and  interesting  applications  concerning  diffusioosmosis 
are  also  available  in  the  literature  [28].  Ionic  solutions  with  a  con¬ 
centration  gradient  of  order  100  kmol  m-4  (=1  M  cm-1 )  along  solid 
surfaces  with  a  zeta  potential  of  order  kT/e  (~25  mV;  e  is  the  charge 
of  a  proton,  k  is  the  Boltzmann  constant,  and  T  is  the  absolute 
temperature)  can  flow  by  diffusioosmosis  at  a  velocity  of  several 
micrometers  per  second. 

The  pores  in  the  membrane  (e.g.,  Nation)  of  a  PEFC  can  be 
simulated  as  straight  nanochannels  with  negatively  charged  walls, 
and  the  hydrogen  ions  generated  at  the  anode  create  a  protonic 
concentration  gradient  in  the  membrane  against  the  cathode;  as 
a  result,  diffusioosmosis  will  contribute  to  the  flow  of  water.  In 


this  work,  we  will  analytically  study  the  effect  of  diffusioosmosis, 
in  comparison  with  that  of  electroosmosis,  on  the  isothermal 
transport  of  water  in  a  fully  hydrated  membrane  whose  shapes  of 
channels  are  maintained  during  the  operation  of  the  fuel  cell.  The 
capillary  channels  are  in  general  of  the  order  nanometer  in  radius, 
and  no  assumption  is  made  about  the  thickness  of  the  electric 
double  layers  adjacent  to  the  capillary  walls  or  the  magnitude  of 
the  zeta  potential  of  the  walls. 


2.  Electrokinetic  flow  in  a  capillary  tube 

In  this  section,  we  consider  the  combined  electroosmotic  and 
diffusioosmotic  flow  of  a  fluid  solution  of  hydrogen  ions  in  a  straight 
capillary  tube  of  radius  R  and  length  L  with  R  <  L  in  the  presence  of 
an  external  electric  field  with  a  constant  magnitude  E  in  the  axial  (z) 
direction,  as  illustrated  in  Fig.  la,  at  the  steady  state.  The  discrete 
nature  of  the  surface  charges,  which  are  negative  and  uniformly 
distributed  over  the  capillary  wall,  is  ignored.  The  prescribed  con¬ 
centration  gradient  Vn°°  of  the  hydrogen  ions  has  a  constant  mag¬ 
nitude  (equal  to  |Vn°°|)  along  the  -z  direction,  where  n°°(z)  is  the 
linear  concentration  distribution  of  the  hydrogen  ions  in  the  bulk 
solution  phase  in  equilibrium  with  the  fluid  inside  the  capillary.  The 
hydrogen  ions  can  diffuse  in  the  capillary,  so  there  exists  no  regular 
osmotic  flow  of  the  solvent.  It  is  assumed  that  L|Vnoo|/noo(0)  <<  1, 
where  z  =  0  is  set  at  the  midpoint  through  the  capillary  axially.  Thus, 
the  variation  of  the  electrostatic  potential  and  ionic  concentration 
in  the  electric  double  layer  (without  co-ions)  adjacent  to  the  capil¬ 
lary  wall  with  the  axial  position  can  be  neglected. 


2.1.  Electrostatic  potential  distribution 

We  first  deal  with  the  electrostatic  potential  distribution  in  the 
ionic  solution  on  a  cross-section  of  the  capillary  tube.  If  xjr{r)  repre¬ 
sents  the  electrostatic  potential  at  a  point  with  distance  r  from  the 
axis  of  the  tube  (where  the  total  electric  potential  equals  xfs-Ez  from 
the  linear  superposition  of  xj/  that  would  exist  in  the  absence  of  the 
imposed  electric  field  and  -Ez  that  arises  from  this  electric  field) 
relative  to  that  in  a  neutral  bulk  solution  and  n(r,  z)  denotes  the 
local  concentration  of  the  hydrogen  ions,  then  Poisson’s  equation 
gives 


l_d  (M\ 

r  dr  y  dr  J 


-~n(r,  0). 
£ 


(1) 


Here  £  is  the  dielectric  permittivity  of  the  ionic  solution,  which  is 
taken  to  be  a  constant.  The  local  concentration  of  hydrogen  ions 
can  also  be  related  to  the  electrostatic  potential  by  the  Boltzmann 
equation: 

n  =  (2) 

where  xj/  =  xj/e/kT  is  the  normalized  electrostatic  potential  distri¬ 
bution. 


Fig.  1.  Geometrical  sketch  for  the  electrokinetic  flow  in  a  straight  nanochannel  due  to  an  applied  electric  field  and  a  prescribed  concentration  gradient  of  hydrogen  ions:  (a) 
a  capillary  tube  and  (b)  a  capillary  slit. 
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Substitution  of  Eq.  (2)  into  Eq.  (1)  leads  to  the 
Poisson-Boltzmann  equation: 


il  (m\ 

r  dr  y  dr  J 


=  -k2  e-* 


(3) 


where  /c  is  the  reciprocal  of  the  Debye  screening  length  defined  by 


k  = 


n°°(0) 


d/2 


skT 


e. 


(4) 


The  boundary  conditions  for  \j/  are 


r  =  0  : 


r  =  R  :  i/s  =  tj, 


(5a) 

(5b) 


where  the  constant  £  =  t;e/kT  is  the  dimensionless  zeta  potential 
at  the  shear  plane  of  the  capillary  wall  adjacent  to  the  fluid  solution 
having  a  uniform  bulk  concentration  n°°( 0)  of  the  hydrogen  ions. 

Introducing  the  variables  y  =  ln(/cr)  and  v=  -^  +  2  ln(/cr),  we 
can  rewrite  Eq.  (3)  as 


d  /dv\2  dey 

dyU y)  d^' 


(6) 


The  integration  of  the  above  equation  with  respect  to  y  subject 
to  boundary  condition  (5a)  leads  to 


2ey +  4. 


(7) 


After  eliminating  the  unreasonable  root,  Eq.  (7)  becomes 


dv 

dy 


=  (2ev  +4)1/2. 


(8) 


By  integrating  the  above  equation  with  respect  to  y  again  subject 
to  boundary  condition  (5b),  the  electrostatic  potential  distribution 
can  be  obtained  as 


\j/  =  2  ln(/cr)  -  In 
where 

7)  =  m  =  2\n 


-2  seek2  [  In 


AT 

‘27=5 


fo 

2 


kR 


+  In 


2  +  a/4  +  2  e-f(/d?)2 


27^2  2  -  a/4  +  2  e-%R)2 


(9a) 


(9b) 


representing  the  dimensionless  electrostatic  potential  on  the  axis 
of  the  capillary  tube. 

If  the  constant  surface  charge  density  cr,  instead  of  the  surface 
potential  £ ,  is  known  at  the  capillary  wall,  the  boundary  condition 
specified  by  Eq.  (5b)  should  be  replaced  by  the  Gauss  condition: 


r/R 

Fig.  2.  Plots  of  the  dimensionless  potential  xj/(r)  in  a  capillary  tube  versus  the  relative 
position  r/R  for  various  values  of  The  solid  and  dashed  curves  denote  the  cases 
kR=  1  and  kR  =  10,  respectively. 


and  kR.  Negative  values  of  £  are  chosen  since  the  capillary  wall  bears 
negative  charges  (e.g.,  resulting  from  sulfonic  acid  (S03_)  groups 
attached  to  the  perfluoro  side  chains  of  a  Nation  carbon-fluorine 
polymer  backbone).  As  expected,  the  value  of  xjr  is  always  nega¬ 
tive  with  its  value  increasing  or  its  magnitude  decreasing  with  a 
decrease  in  r/R  when  the  electric  double  layer  adjacent  to  the  cap¬ 
illary  wall  is  relatively  thick  (say,  with  /eft  <  1 )  for  a  constant  value  of 
£  at  the  wall.  On  the  other  hand,  when  the  double  layer  is  relatively 
thin  (say,  with  /eft  >2),  the  value  of  \j/  increases  with  a  decrease  in 
r/ft  from  the  given  dimensionless  zeta  potential  £  at  the  wall  (with 
r/ft  =  1 )  and  becomes  positive  after  r/ft  is  smaller  than  certain  val¬ 
ues  [29],  since  the  fluid  around  the  axis  of  the  capillary  tube  still 
contains  hydrogen  ions  only  and  the  negatively  charged  wall  has 
little  influence  on  the  fluid  in  this  region.  For  a  specified  value  of 
/eft,  the  value  of  ^  at  any  position  r/ft  decreases  with  a  decrease 
in  the  value  of  £.  It  can  be  seen  that  the  value  of  \j/  -  £  is  always 
positive  and  increases  with  an  increase  in  the  magnitude  of  £  for 
given  values  of  r/ft  and  /eft.  The  increase  of  xj/  -  £  with  decreasing  r/ft 
(which  is  relatively  sharp  near  the  capillary  wall)  becomes  gentler 
when  either  /eft  or  the  magnitude  of  £  is  smaller. 

2.2.  Fluid  velocity  distribution 


dxjr  ea 
dr  -  ekT' 


(10) 


The  solution  for  \fr  given  by  Eqs.  (9a)  and  (9b)  still  holds  for  this 
condition,  with  the  relation  between  £  and  cr  for  an  arbitrary  value 
of  /eft  as 


?  =  ln 


2  (/d?)2 

{{eoR/ekT)  -  2)2  -  4 


(11) 


We  now  consider  the  steady  electrokinetic  flow  of  the  hydrogen 
ionic  solution  in  the  capillary  tube.  The  momentum  balances  on  the 
Newtonian  fluid  in  the  r  and  z  directions  give 


dp  d^/  n 

-f+en-f-  =  0, 
dr  dr 


(12a) 


5—  (r  —  ) 1 
r  dr  V  dr  / 


dp 

dz 


-  enE, 


(12b) 


For  the  hydrogen  ionic  solution  in  a  capillary  tube  with  a  given 
radius,  the  above  equation  predicts  that  the  magnitude  of  f 
decreases  with  an  increase  in  /eft  for  the  case  of  constant  surface 
charge  density  and  the  magnitude  of  cr  increases  with  an  increase 
in  /eft  for  the  case  of  constant  surface  potential. 

Fig.  2  shows  the  result  of  the  dimensionless  electrostatic  poten¬ 
tial  \j/  calculated  from  Eq.  (9a)  and  (9b)  as  functions  of  the  relative 
position  r/ft  for  various  values  of  the  dimensionless  parameters  £ 


where  u(r)  is  the  fluid  velocity  in  the  positive  z  direction  (axial 
direction  of  decreasing  ionic  concentration),  p(r,  z)  is  the  pressure 
distribution,  r\  is  the  viscosity  of  the  fluid,  and  E  is  the  imposed 
electric  field  in  the  z  direction.  The  boundary  conditions  for  u  are 


r  =  0  : 


r  =  ft  : 


?=0, 

dr 


(13a) 

(13b) 


u  =  0. 
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r/R 


Fig.  3.  Plots  of  the  dimensionless  pressure  (p  —po)ln°°kT  in  a  capillary  tube  versus 
the  relative  position  r/R  for  various  values  of  £.  The  solid  and  dashed  curves  denote 
the  cases  kR  =  1  and  kR  =  10,  respectively. 


After  the  substitution  of  Eq.  (2)  into  Eq.  (12a),  the  pressure  dis¬ 
tribution  can  be  determined  as 

p  =  Po  +  n°°(z)fcr[e-^  -  e-'H  (14) 

where  p0  is  the  pressure  on  the  axis  of  the  capillary  tube,  which 
is  a  constant  in  the  absence  of  applied  pressure  gradients  and  the 
electrostatic  potential  distribution  \js(r)  is  given  by  Eqs.  (9a)  and 
(9b).  Note  that  the  pressure  distribution  in  the  r  and  z  directions 
results  from  the  variation  of  the  electrostatic  potential  ^  with  r 
and  the  variation  of  the  bulk  concentration  n°°  with  z,  respectively. 

Fig.  3  illustrates  the  result  of  the  normalized  pressure 
C p-Po)ln°°kT obtained  from  Eqs.  (14),  (9a)  and  (9b)  as  functions  of 
the  relative  coordinate  r/R  for  several  values  of  the  dimensionless 
parameters  £  and  kR.  It  can  be  seen  that  this  normalized  pressure 
is  always  positive  and  increases  with  an  increase  in  r/R  irrespective 
of  the  values  of  f  and  kR.  This  behavior  is  understood  by  an  exam¬ 
ination  of  Eq.  (12a)  for  the  fluid  momentum  balance  in  the  radial 
direction  and  knowing  that  the  value  of  the  electrostatic  poten¬ 
tial  \j/  (which  is  negative  at  the  capillary  wall)  decreases  with  an 
increase  in  r/R.  Accordingly,  the  value  of  (p  -p0)ln°°kT  increases 
with  an  increase  in  the  magnitude  of  £  for  given  values  of  r/R  and 
kR,  and  its  increase  with  r/R  (which  is  relatively  sharp  near  the  cap¬ 
illary  wall)  becomes  gentler  when  either  kR  or  the  magnitude  of  £ 
is  smaller.  Eq.  (14)  also  predicts  that  the  axial  gradient  of  the  excess 
pressure  in  the  electric  double  layer  is  in  the  same  direction  as  that 
of  the  ionic  concentration  gradient. 

Substituting  the  ionic  concentration  distribution  of  Eq.  (2)  and 
the  pressure  profile  of  Eq.  (14)  into  Eq.  (12b),  and  then  perform¬ 
ing  the  integration  with  respect  to  r  twice  subject  to  the  boundary 
conditions  in  Eqs.  (13a)  and  (13b),  we  obtain 

/r/R  n  r/R 

-  J  J[(e-*>-e-*)N-e-*E ]dr-dr~,  (15) 

where  the  electrostatic  potential  distribution  \jr  is  given  by  Eqs.  (9a) 
and  (9b), 

(i6> 

tjk 

which  is  a  characteristic  velocity  of  the  electrokinetic  flow  of 
the  fluid,  E  =  Ee/KkT  represents  a  dimensionless  strength  of 
the  applied  electric  field  to  produce  the  electroosmosis,  and 


N  =  |Vnoo|//cnoo(0)  denotes  a  dimensionless  strength  of  the  ionic 
concentration  gradient  to  cause  diffusioosmosis  (more  accu¬ 
rately,  chemiosmosis  only).  For  some  typical  values,  E &  250  V  m-1 , 
n°7|Vn°°|  «  1(T4  m,  and  /c-1  ^10-8  m,  then  both  E  and  N  have  the 
order  10-4.  Eq.  (15)  indicates  that  the  total  fluid  velocity  is  a  super¬ 
position  of  the  electroosmotic  and  diffusioosmotic  velocities  and, 
under  the  condition  of  N  =  E,  the  effect  of  diffusioosmosis  is  weaker 
than  that  of  electroosmosis  by  a  factor  1  - 

The  average  fluid  velocity  (u)  over  a  cross-section  of  the  capillary 
tube  is  defined  as 

,",-2/o,“(054 

which  can  be  obtained  by  numerical  integration  after  the  substitu¬ 
tion  of  Eq.  (15). 

2.3.  Results  of  fluid  velocity  and  discussion 

The  electroosmotic  and  diffusioosmotic  velocities  of  the  hydro¬ 
gen  ionic  solution  in  the  capillary  tube  are  given  by  Eq.  (15),  in 
which  the  dual  integration  can  be  performed  numerically,  setting 
N  =  0  and  E  =  0,  respectively.  These  velocities  in  nondimensional 
mobility  form  are  plotted  versus  the  normalized  coordinate  r/R  in 
Figs.  4  and  5  for  several  values  of  the  parameters  f  and  kR.  It  can 
be  seen  that  both  the  dimensionless  electroosmotic  mobility  u/U*E 
and  diffusioosmotic  mobility  u/U*N  decrease  monotonically  with 
an  increase  in  r/R  from  maxima  at  the  axis  of  the  capillary  ( r/R  =  0) 
to  zero  at  the  capillary  wall  (r/R  =  1)  due  to  the  no-slip  condition 
at  the  wall.  The  fluid  flows  in  the  direction  of  the  electric  field  or 
decreasing  ionic  concentration  (u  is  positive).  Both  the  electroos¬ 
motic  and  diffusioosmotic  mobilities  are  monotonically  increasing 
functions  of  kR  and  the  magnitude  of  £.  It  implies  that  a  modifi¬ 
cation  of  the  membrane  material  in  a  fuel  cell  with  a  smaller  pore 
radius  or  a  lower  magnitude  of  zeta  potential  could  result  in  less 
water  transport  and  easier  maintenance  of  the  fully  hydrated  envi¬ 
ronment.  Since  the  aging  of  the  membrane  and  the  deformation  of 
its  nanopores  during  fuel  cell  operations  often  change  the  values 
of  kR  and  £,  these  factors  could  affect  the  consistency  of  the  water 
transport  as  well. 


Fig.  4.  Plots  of  the  dimensionless  electroosmotic  mobility  u/U*E  in  a  capillary  tube 
as  calculated  from  Eq.  (15)  with  N  =  0  versus  the  relative  position  r/R  for  vari¬ 
ous  values  of  The  solid  and  dashed  curves  denote  the  cases  kR  =  1  and  kR  =  2, 
respectively. 
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Fig.  5.  Plots  of  the  dimensionless  diffusioosmotic  mobility  u/U*E  in  a  capillary  tube 
as  calculated  from  Eq.  (15)  with  E  =  0  versus  the  relative  position  r/R  for  various  val¬ 
ues  of  £.  The  solid  and  dashed  curves  denote  the  cases  kR  =  1  and  kR  =  2,  respectively. 

In  Fig.  6,  the  dimensionless  average  electrokinetic  velocity 
(■ u)/U*E  of  the  ionic  solution  in  the  capillary  tube  calculated  using 
Eq.  (17)  is  plotted  as  functions  of  the  parameter  kR  for  some  values 
of  the  parameters  N/E  and  f.  The  contribution  from  diffusioosmo- 
sis  to  the  total  fluid  flow  increases  with  the  ratio  N/E,  as  expected, 
and  is  quite  significant  (more  than  20%)  as  this  ratio  is  greater  than 
1  / 2.  In  accordance  to  the  electroosmotic  and  diffusioosmotic  veloc¬ 
ity  profiles  illustrated  in  Figs.  4  and  5,  the  average  fluid  velocity 
increases  with  an  increase  in  kR  and  the  magnitude  of  f  at  any 
value  of  N/E. 


two  identical  parallel  plates  with  a  separation  distance  2 h  in  the 
presence  of  a  uniform  external  electric  field  E  in  the  tangential 
z  direction,  as  shown  in  Fig.  lb.  The  imposed  concentration  gra¬ 
dient  Vn°°  of  the  hydrogen  ions  is  along  the  -z  direction.  Owing 
to  the  planar  symmetry  of  the  system,  we  need  consider  only  the 
half  region  0  <y  <  h,  where  y  is  the  distance  from  the  median  plane 
between  the  slit  walls  in  the  normal  direction.  The  analysis  for  this 
case  is  similar  to  that  for  the  case  of  electrokinetic  flow  in  a  capillary 
tube  presented  in  the  previous  section. 


3.1  Electric  potential  distribution 


For  the  ionic  solution  in  the  capillary  slit,  let  ^(y)  be  the 
electrostatic  potential  at  the  position  y  and  n(y,  z)  be  the  local  con¬ 
centration  of  the  hydrogen  ions.  Flere,  Poisson’s  equation  becomes 


d2\[/ 

dy2 


-e-n(y,  0), 


(18) 


where  z  =  0  is  set  at  the  midpoint  through  the  capillary.  Substitution 
of  Eq.  (2)  into  the  above  equation  results  in 


d  2\j/ 
dy2 


-k2 


(19) 


where  the  Debye  screening  parameter  k  is  defined  by  Eq.  (4).  For 
the  case  of  constant  surface  potential  £,  the  boundary  conditions 
for  are 


y  =  0  : 


y  =  h  :  &  = 


(20a) 

(20b) 


The  solution  to  Eq.  (19)  subject  to  Eqs.  (20a)  and  (20b)  can  be 
obtained  as 


\js  =  \j/Q  +  2  lncos[C  -  (/ch  -  /cy)B], 


(21) 


3.  Electrokinetic  flow  in  a  capillary  slit 

In  this  section,  we  examine  the  steady  electroosmosis  and  diffu- 
sioosmosis  of  a  hydrogen  ionic  solution  in  a  nanochannel  between 


Fig.  6.  Plots  of  the  dimensionless  average  electrokinetic  velocity  (u)  /U*E  in  a  cap¬ 
illary  tube  as  calculated  from  Eq.  (17)  versus  the  electrokinetic  radius  kR  for  various 
values  of  N/E.  The  solid  and  dashed  curves  denote  the  cases  £  =  -1  and  £  =  -2, 
respectively. 


where 

B  =  e-«A o/2/V2,  (22a) 

C  =  tan”1  (22b) 


and  now  =  \j/{ 0)  =  f  -  2  In  cos  C,  representing  the  potential  at 
the  median  plane  between  the  slit  walls. 

For  the  case  of  constant  surface  charge  density  o ,  the  boundary 
condition  given  by  Eq.  (20b)  is  replaced  by 


y  =  h: 


d \j/  ecr 
dy  -  £kT  ’ 


(23) 


and  the  solution  for  \jr  given  by  Eq.  (21 )  is  still  valid,  with 

1  = 1  +  i^rcoth('ch5-  (24) 

For  the  ionic  solution  in  a  capillary  slit  with  a  given  separation 
distance  between  its  walls,  the  magnitude  of  £  is  a  decreasing  func¬ 
tion  of  Kh  for  a  constant  value  of  o  and  the  magnitude  of  a  is  an 
increasing  function  of  /ch  for  a  constant  value  of  if. 

The  numerical  values  of  the  dimensionless  potential  \j/  cal¬ 
culated  from  Eq.  (21)  are  plotted  in  Fig.  7  as  functions  of  the 
dimensionless  coordinate  y/h  for  several  values  of  the  dimension¬ 
less  parameters  £  and  \ch.  Similar  to  the  case  of  a  capillary  tube 
shown  in  Fig.  2,  the  value  of  \jr  is  always  negative  with  its  magni¬ 
tude  decreasing  with  a  decrease  in  y/h  when  the  electric  double 
layers  adjacent  to  the  walls  are  relatively  thick  for  a  constant  mag¬ 
nitude  of  and  it  increases  with  a  decrease  in  y/h  from  the  given  l 
at  the  walls  and  becomes  positive  after  y/h  is  smaller  than  certain 
values  when  the  double  layers  are  relatively  thin.  The  value  of  \j/  at 
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Fig.  7.  Plots  of  the  dimensionless  potential  \js(y)  in  a  capillary  slit  versus  the  relative 
position  y/h  for  various  values  of  £ .  The  solid  and  dashed  curves  denote  the  cases 
Kh  =  1  and  ich  =  10,  respectively. 


any  position  decreases  with  a  decrease  in  the  values  of  £  and  Kh. 
Since  the  adjacent  wall  area  per  unit  volume  of  the  fluid  is  smaller 
in  a  slit  than  in  a  tube,  the  decay  of  the  electrostatic  potential  from 
the  capillary  walls  is  faster  in  a  slit  than  in  a  tube. 


3.2.  Fluid  velocity  distribution 


The  momentum  equations  for  the  steady  flow  of  the  solution  of 
hydrogen  ions  in  a  capillary  slit  with  an  applied  electric  field  and 
ionic  concentration  gradient  parallel  to  the  slit  walls  are 


dp  d^  _ 

/+en/  =0, 
dy  d y 

d  2u  dp 
r\~r^  =  q-  ~enE, 
dy2  dz 


(25a) 

(25b) 


where  u(y)  is  the  fluid  velocity  profile  in  the  tangential  direction  of 
decreasing  ionic  concentration,  p(y,  z)  is  the  pressure  distribution, 
and  E  is  the  applied  electric  field.  Here,  the  boundary  conditions  for 
u  are 

y  =  o  :  ^  =  0,  (26a) 

y  =  h  :  u  =  0.  (26b) 


After  the  substitution  of  Eq.  (2)  into  Eq.  (25a),  the  pressure  distri¬ 
bution  can  still  be  determined  as  the  form  of  Eq.  ( 14),  in  which  p0  is 
the  pressure  on  the  midplane  between  the  slit  walls  and  the  electro¬ 
static  potential  distribution  ijs(y)  is  given  by  Eq.  (21 ).  Fig.  8  illustrates 
the  result  of  the  normalized  pressure  (p  -p0)ln°°kT  obtained  from 
Eqs.  (14)  and  (21 )  as  functions  of  the  relative  coordinate  y/h  for  var¬ 
ious  values  of  the  dimensionless  parameters  £  and  Kh.  Analogous 
to  the  result  in  the  case  of  a  capillary  tube  in  Fig.  3,  (p  -po)ln°°kT  is 
positive  and  an  increasing  function  of  y/h  for  any  specified  values 
of  £  and  Kh.  This  normalized  pressure  increases  with  an  increase  in 
the  magnitude  of  £  for  given  values  of  y/h  and  Kh,  and  its  increase 
with  y/h  becomes  gentler  when  either  /ch  or  the  magnitude  of  If 
is  smaller.  Since  the  decay  of  the  electrostatic  potential  from  the 


Fig.  8.  Plots  of  the  dimensionless  pressure  (p -p0)ln°°kT  in  a  capillary  slit  versus 
the  relative  position  y/h  for  various  values  of  £.  The  solid  and  dashed  curves  denote 
the  cases  Kh  =  1  and  Kh  =  10,  respectively. 


capillary  walls  is  faster  in  a  slit  than  in  a  tube,  the  variation  of  the 
pressure  with  the  lateral  position  is  stronger  in  a  slit  than  in  a  tube, 
as  Eq.  (25a)  predicts. 

The  fluid  velocity  u  can  be  solved  from  Eqs.  (25b),  (26a)  and 
(26b)  by  the  substitution  of  the  ionic  concentration  distribution  of 
Eq.  (2)  and  the  pressure  profile  of  Eq.  (14),  with  a  form  similar  to 
Eq.  (15): 

11  fy/h  fy/h  -  -  -  V  V 

(L={lchfj  J  Ue-+°-e~+)N-e-+E]dldl,  (27) 

where  U*  is  defined  by  Eq.  (16).  After  substituting  the  electrostatic 
potential  distribution  of  Eq.  (21 )  and  integrating  with  respect  to  y 


y/h 

Fig.  9.  Plots  of  the  dimensionless  electroosmotic  mobility  u/U*E  in  a  capillary  slit  as 
calculated  from  Eq.  (28)  with  N  =  0  versus  the  relative  positiony/h  for  various  values 
of  £.  The  solid  and  dashed  curves  denote  the  cases  Kh  =  1  and  Kh  =  2,  respectively. 
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Fig.  10.  Plots  of  the  dimensionless  diffusioosmotic  mobility  u/U*N  in  a  capillary 
slit  as  calculated  from  Eq.  (28)  with  E  =  0  versus  the  relative  position  y/h  for  var¬ 
ious  values  of  £.  The  solid  and  dashed  curves  denote  the  cases  ich  =  1  and  Kh  =  2, 
respectively. 

twice,  the  above  equation  becomes 
£  =  -m)2-(Ky)2]B2N 

+  2  {lnC°S[C-c(0fc^y)B|  -  (*/.  -  Ky)Btan(C  -  chB)} 

x  (£  +  N),  (28) 

where  the  dimensionless  parameters  B  and  C  are  defined  by  Eqs. 
(22a)  and  (22b). 

The  average  fluid  velocity  (u)  over  a  cross-section  of  the  slit  can 
be  expressed  as 


which  can  be  obtained  by  numerical  integration  after  the  substitu¬ 
tion  of  Eq.  (28). 


3.3.  Results  of  fluid  velocity  and  discussion 

The  electroosmotic  velocity  of  the  hydrogen  ionic  solution  in  the 
capillary  slit  is  given  by  Eq.  (28)  setting  N  =  0,  whereas  the  diffu¬ 
sioosmotic  velocity  can  be  obtained  by  the  same  equation  taking 
E  =  0.  In  Figs.  9  and  10,  these  electrokinetic  velocities  in  dimension¬ 
less  mobility  form  are  plotted  versus  the  normalized  coordinate 
y/h  for  several  values  of  the  parameters  f  and  /ch.  Analogous  to  the 
case  of  a  capillary  tube  discussed  in  the  previous  section,  both  the 
dimensionless  electroosmotic  mobility  u/U*E  and  diffusioosmotic 
mobility  u/U*N  decrease  monotonically  with  an  increase  in  y/h  and 
vanish  at  the  no-slip  capillary  walls.  The  fluid  flows  toward  the  end 
with  lower  electric  potential  or  ionic  concentration.  Both  of  these 
electrokinetic  velocities  are  monotonically  increasing  functions  of 
Kh  and  the  magnitude  of  f. 

The  dimensionless  average  electrokinetic  velocity  (u) /U*E  of  the 
ionic  solution  in  the  capillary  slit  calculated  using  Eq.  (29)  is  plot¬ 
ted  in  Fig.  11  as  functions  of  /ch  for  some  values  of  the  parameters  \ 
and  N/E.  Again,  the  contribution  from  diffusioosmosis  to  the  total 
fluid  flow  increases  with  the  ratio  N/E  and  is  quite  significant  as 
this  ratio  is  greater  than  1/2.  In  accordance  to  the  electroosmotic 


and  diffusioosmotic  velocity  profiles  illustrated  in  Figs.  9  and  10, 
the  average  fluid  velocity  increases  with  an  increase  in  /ch  and  the 
magnitude  of  f  at  any  value  of  N/E.  Since  the  viscous  retardation 
exerted  on  the  fluid  flow  by  the  capillary  walls  decreases  with  a 
decrease  in  the  adjacent  wall  area  per  unit  volume  of  the  fluid,  the 
normalized  magnitudes  of  the  electroosmotic  velocity,  diffusioos¬ 
motic  velocity,  and  total  electrokinetic  velocity  in  a  slit  is  greater 
than  their  corresponding  magnitudes  in  a  tube. 

4.  Concluding  remarks 

An  analytical  study  of  the  effect  of  diffusioosmosis,  in  compar¬ 
ison  with  that  of  electroosmosis,  on  the  transport  of  water  in  the 
fully  hydrated  membrane  of  a  PEFC  at  the  steady  state  is  presented 
in  this  work.  The  reactant  gases  are  fed  at  the  same  pressure  so 
that  no  external  pressure  gradient  exists  in  the  membrane  and  the 
excess  water  produced  by  the  chemical  reaction  at  the  cathode  is 
drained.  A  capillary  tube  or  slit  with  negatively  charged  walls  is 
chosen  to  model  the  interstitial  nanochannels  in  the  membrane.  It 
is  assumed  that  the  fluid  is  only  slightly  nonuniform  in  the  concen¬ 
tration  of  hydrogen  ions  along  the  axial  direction  of  the  capillary, 
but  no  assumption  is  made  about  the  thickness  of  the  electric  dou¬ 
ble  layers  adjacent  to  the  capillary  walls.  The  capillary  walls  may 
have  either  a  constant  surface  potential  or  a  constant  surface  charge 
density  of  an  arbitrary  quantity.  By  solving  the  Poisson-Boltzmann 
equation  and  the  modified  Navier-Stokes  equation  applicable  to 
the  system,  the  electrostatic  potential  distribution  and  the  pressure 
profile  under  the  influence  of  a  constant  prescribed  concentration 
gradient  of  hydrogen  ions  are  determined  analytically.  The  results 
for  the  local  and  averaged  electrokinetic  velocities  of  the  fluid  in  the 
capillary  are  presented  in  detail,  and  the  effect  of  diffusioosmosis 
on  the  water  transport  in  the  membrane  of  a  PEFC  is  found  to  be 
significant,  especially  when  the  macroscopic  electric  field  is  low. 

Although  we  discuss  the  diffusioosmotic  effect  in  the  previ¬ 
ous  sections  mainly  on  PEFCs,  with  protons  as  the  moving  ions, 
the  same  analysis  can  also  be  applied  for  other  fuel  cells  with  a 
similar  mechanism.  For  examples,  in  solid-state  alkaline  fuel  cells, 
hydroxide  ions  migrate  in  anion-exchange  membranes,  whereas  in 
molten  carbonate  fuel  cells,  carbonate  ions  move  in  chemically  inert 


Fig.  11.  Plots  of  the  dimensionless  average  electrokinetic  velocity  ( u)/U*E  in  a  cap¬ 
illary  slit  as  calculated  from  Eq.  (29)  versus  the  electrokinetic  radius  Kh  for  various 
values  of  N/E.  The  solid  and  dashed  curves  denote  the  cases  £  =  -1  and  £  =  -2, 
respectively. 
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ceramic  matrices.  The  major  difference  is  that  the  positive  hydro¬ 
gen  ions  in  the  membranes  are  replaced  by  some  negative  ions,  and 
the  direction  of  the  electrokinetic  flow  of  the  fluid  will  be  reversed. 
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